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Abstract 

Using the helicity method we derive complete formulas for the joint angular decay dis- 
tributions occurring in semileptonic hyperon decays including lepton mass and polarization 
effects. Compared to the traditional covariant calculation the helicity method allows one to 
organize the calculation of the angular decay distributions in a very compact and efficient 
way. In the helicity method the angular analysis is of cascade type, i.e. each decay in the 
decay chain is analyzed in the respective rest system of that particle. Such an approach is 
ideally suited as input for a Monte Carlo event generation program. As a specific example 
we take the decay 2° — > S + + l~ + v\ (l~ = e~,/i~) followed by the nonleptonic decay 
S + — > p + 7T° for which we show a few examples of decay distributions which are generated 
from a Monte Carlo program based on the formulas presented in this paper. All the results 
of this paper are also applicable to the semileptonic and nonleptonic decays of ground state 
charm and bottom baryons, and to the decays of the top quark. 



1 Introduction 



Semileptonic hyperon decays have traditionally been analyzed in the rest frame of the 
decaying parent hyperon using fully covariant methods based on either four-component 
Dirac spinor methods [HEHSlHJElElEjoron two-component Pauli spinor methods [8|l9l[T0]. 
The latter method is particularly well suited for an implementation of a zero- [8] or a 
near zero-recoil [9} \TU\ approximation. In the present paper we employ helicity methods 
to analyze semileptonic hyperon decays. In the muonic mode it is quite important to 
incorporate lepton mass effects in the analysis since e.g. in the decay 2° — > S + + l~ + v\ 
the mass difference between the parent and daughter hyperon M=o — = (1314.83 — 
1189.37) MeV = 125.46 MeV is comparable to the muon mass m M = 105.658 MeV. The 
analysis proceeds in cascade fashion where every decay in the decay chain is analyzed in its 
respective rest frame. For the semileptonic decays 5° — > J 9 + 7r °) + W / fr-shcii(^ ^ 

(I = e~, treated in this paper this means that the decay 2° — > £ + + W~ s _ shell is analyzed 
in the S° rest frame whereas the decays S + — ► p + n° and W~ s _ shcll — > l~ + vi are analyzed 
in the S + and W~ s _ shell rest frames, respectively. In this way one obtains exact decay 
formulas with no approximations which are quite compact since they can be written in a 
quasi-factorized form. 

Cascade-type analysis were quite popular some time ago in the strong interaction sector 
when analyzing the decay chains of the strong interaction baryonic and mesonic resonances 
(see e.g. [ID [131 III]). I n the weak interaction sector cascade-type analysis were applied 
before to nonleptonic decays [H2, [16j [TFJ [181 EEl EH EH, [22], to semileptonic decays of heavy 
mesons and baryons [HI HH HH [23l [2H [25l [261 12ZI [281 ISH [30] , and to rare decays of heavy 
mesons [301 131] and heavy baryons [32]. A new feature appears in semileptonic decays 
compared to nonleptonic decays when one includes lepton mass effects. In this case one 
has new interference contributions coming from the time-components of the vector and 
axial vector currents interfering with the usual three-vector components of the currents 
(see e.g. [HIES]). 

The results for the angular decay distributions in the semileptonic decays of heavy 
baryons given e.g. in [18], [28] can in fact be directly transcribed to the hyperon sector q 
However, the presentation in [TH [28] is rather concise and concentrates on results for 
angular decay distributions and their analysis rather than presenting the details of their 
derivations. In order to make the results more reproducible we decided to include the details 
of the derivations in this paper. This will enable the interested reader to e.g. convert the 
results of [18], which were derived for the (Z + , vi) case, to the v{) case discussed in this 

1 Our approach is similar in spirit to the approach of [33j which also used a cascade-type helicity analysis 
to analyze the semileptonic decay of a polarized hyperon using the lepton side as a polarization analyzer. 
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paper, or, to derive angular decay formula involving the decay E + — > p + 7 instead of the 
decay S + — ► p + tt° treated in this paper. At the same time we decided to recalculate all 
relevant decay formulas in order to provide another independent check of their correctness. 
In this way we discovered one error in [28] and two errors in [18J. 

In the simulation of semileptonic hyperon decays including the /i-mode it is important 
to have a reliable and tested Monte Carlo (MC) program. Since hyperons are produced 
with nonzero polarization the MC program should also include polarization effects of the 
decaying parent hyperon. One of the motivations of starting this project was the fact 
that such a general purpose MC event generator has not been available up to now. Such 
an event generator should prove to be quite useful in the analysis of the huge amount of 
data on semileptonic hyperon decays that has been collected by the KTeV and the NA48 
collaborations. We wrote and tested such a MC event generator based on the formulas 
written down in this paper. The present paper can be viewed as a documentation of 
the theoretical spin-kinematical input that goes into the MC program and, for the sake 
of reproducibility, the paper also describes how to derive the angular decay distributions 
entering the MC. 

Although we frequently refer to the specific semileptonic cascade decay S° — > 
p + it ) + W^_ shell (— > l~ + 9i) the spin-kinematical analysis presented in this paper is 
quite general and can be equally well applied to the semileptonic decays of heavy charm 
and bottom baryons, and for that matter, also to the semileptonic decay of the top quark. 
In order to facilitate such further applications we have always included the necessary sign 
changes when going from the (l~ , vj) case to the (7 + , vj) case as occurs e.g. in the semilep- 
tonic hyperon decay S + — > A + e + + v e , in semileptonic c — > s charm baryon decays or in 
semileptonic top quark decays pU [35], [36j [37] . When sign changes are indicated the upper 
sign will always refer to the (l~ , vj) case, which is the main concern of this paper, whereas 
the lower sign will refer to the (7 + , v\) case. We also mention that we have always assumed 
that the amplitudes are relatively real and have therefore dropped azimuthal correlation 
contributions coming from the imaginary parts. Put in a different language this means 
that we have not considered T-odd contributions in our angular analysis which could re- 
sult from final state interaction effects or from truly CP-violating effects. By keeping the 
imaginary parts in the azimuthal correlation terms one can easily write down the relevant 
T-odd contributions if needed by using the formulas of this paper. This is discussed for a 
specific example in Appendix ICl 

The paper is structured as follows. In Sec. [2] we introduce the helicity amplitudes and 
relate them to a standard set of invariant form factors. In order to estimate the size of 
the helicity amplitudes for the H° — > S + current-induced trasition we provide some simple 
estimates for the invariant form factors and their g 2 -dependence which we shall refer to 
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as the minimal form factor model. In Sec. [3] we derive the unpolarized decay rate written 
in terms of bilinear forms of the helicity amplitudes. Sec. H] contains some numerical 
results on branching ratios, the rate ratio T(e)/T(/j,) and a lepton-side forward-backward 
asymmetry. In Sec. [5] we discuss single spin polarization effects including spin-momentum 
correlation effects between the polarization of the parent baryon and the momenta of the 
decay products. Sec. [6] treats momentum-momentum correlations between the momenta 
of the decay products in the cascade decay S° — > > p + 7r°) + W^_ shell (— > l~ + v{) for 
an unpolarized 5°. In Sec. [7] we present a few sample distributions generated from the MC 
program written by us. Sec. M contains our summary and our conclusions. 

We have collected some technical material in the appendices. In Appendix [A] we re- 
count how the two-body decay of a polarized particle is treated in the helicity formalism. 
This two-body decay enters as a basic building block in our quasi-factorized master formu- 
lae in the main text which describe the various cascade-type angular decay distributions 
presented in this paper. In Appendix [B] we list explicit forms of the Wigner's <i J -function 
for J = 1/2 and J = 1, which are needed in the present application. In Appendix Owe go 
through a specific example and identify a specific T-odd term in the joint angular decay 
distribution written down in Sec. The example is easily generalized to other cases. In 
Appendix [D] we finally list the full five-fold angular decay distribution for the cascade de- 
cay S° -> £+(-> P + tt°) + W- ff _ 

-shell (~~ * ^ + u i) f° r a polarized parent hyperon H°. The full 
five-fold angular decay distribution reduces to the decay distributions listed in the main 
text after integration or after setting the relevant parameters to zero. 



2 The helicity amplitudes 



The momenta and masses in the semileptonic hyperon decays are denoted by -BiQoi, Mi) — > 
£> 2 (jo 2 ,M 2 ) + l(pi,mi) + ui(p u ,0). For the hadronic transitions described by the helicity 
amplitudes it is not necessary to distinguish between the cases and (7 + ,z/;). The 

matrix elements of the vector and axial vector currents JY' between the spin 1/2 states 
are written as 



M, 



Mi 



v 



{B 2 \J v \B x )=u 2 (p 2 ) 



{B 2 \J*\B X ) = u 2 (p 2 ) 



F 2 A (q 2 
Mi 



EM 

Mi 
Mi 



ui{pi) , 
75«i0i) , 



(1) 
(2) 



where q = p\ — p 2 is the four-momentum transfer. As in [7] we take a^ u = \{^^ u — JuJn) 
and 75 = — ( 1 J . The other 7 matrices are defined as in Bjorken-Drell. 
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Next we express the vector and axial vector helicity amplitudes H^ A w (A 2 = ±1/2; 
Xw = t, ±1, 0) in terms of the invariant form factors, where and A 2 are the helicity 
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components of the W fr- S heii and the daughter baryon, respectively. Since lepton mass 
effects are taken into account in this paper we need to retain the time-component "t" of 
the four-currents J^' A - Concerning the transformation properties of the four components 
of the currents one notes that, in the rest frame of the VFoff-sheii (<f — 0), the three space- 
components Xw = ±1,0 transform as J = 1 whereas the time-component transforms as 
J = 0. In this paper we use a short-hand notation Xw = t, ±1,0 for Xw = 0(J = 
0); ±1, 0(J = 1). Whenever we write Xw = t this has to be understood as Xw = 0(J = 0). 
One then needs to calculate the expressions 

H H W = Ml>\X 2 )e^{X w ) . (3) 

We do not explicitly denote the helicity (m-quantum number) of the parent hyperon Ai in 
the helicity amplitudes since Ai is fixed by the relation Ai = A2 — Xw- It is very important 
to detail the phase conventions when evaluating the expression in Eq. ([3]). This is because 
the angular decay distributions to be discussed later on contain interference contributions 
between different helicity amplitudes which depend on the relative signs of the helicity 
amplitudes. We shall work in the rest frame of the parent baryon B\ with the daughter 
baryon B2 moving in the positive z-direction. The baryon spinors are then given by [38] 

«i(±±,j>i) = V2mJ *± ) , (4) 
where x+ — [ n ) an< ^ X- = \ ? J are the usual Pauli two-spinors. For the four polar- 



j ^ 1 / 

ization four-vectors of the currents we have [38 



^(t) = -^=(g o ;0,0,-p) 



6^(±1) = -L(0;±l,-i,0) 



^(0) = (p; 0,0, -q ) , (5) 



where g M = (go; 0, 0, —p) is the momentum four-vector of the off-shell gauge boson PVofr-shcii 
in the B\ rest frame. The energy of the off-shell VT-boson go and the magnitude of three- 
momentum of the daughter baryon B 2 (or the off-shell VT-boson) p in the rest system of 
the parent baryon B\ are given by 

* = ^-(M 1 2 -M| + g 2 ), (6) 
P = IP2I = i^VQ^ , (7) 
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where 

Q± = (M 1 ±M 2 ) 2 -q 2 . (8) 

The bar over the polarization four-vectors reminds one that the m quantum numbers of 
the currents are quantized along the negative z-direction. They are obtained from the 
polarization four-vectors quantized along the positive z-axis by a 180° rotation around the 
y-axis (see [38]). Using the spinors in Eq. (j4]) and the polarization vectors Eq. ([5]) one 
obtains following vector helicity amplitudes (Ai = A 2 — \w) 

HY_ X = ^2QI(-F^-{M 1 +M 2 )/M 1 fA , 

HY = ^((M 1 + M 2 )if + q 2 /M 1 F 2 v y (9) 
and axial vector helicity amplitudes 

Hf t = ^(^-(M 1 + M 2 )F 1 A + q 2 /M 1 F 3 A y 

Hf = ^^-( Ml ~M 2 )F 1 A + q 2 /M 1 F 2 A y (10) 



From parity or from an explicit calculation one has 

H v - H v 

— A2,— \w A2,A^ 

H—\2,~\w H\2,Xw (H) 

When discussing the semileptonic transitions close to the zero recoil point it is ad- 
vantageous to make use of the velocity transfer variable uo = v\ ■ v 2 where t> M = p^/M. 
The velocity transfer variable can be expressed in terms of the usual momentum transfer 
variable q 2 . One has 

u = {Ml + Ml - q 2 )/(2M l M 2 ) , (12) 

or 

q 2 = (Mi -M 2 f -2M 1 M 2 (uj- 1) . (13) 

The maximal and minimal values of the velocity transfer variable u are uj ma , x = (Mf + M% — 
mf)/(2MiM 2 ) and u> m i n = 1, respectively. The minimal value u m i n = 1 is referred to as the 
zero recoil point since this is the point where the recoiling daughter baryon has no three- 
momentum. For the variables Q± defined in Sec. [2] one finds Q± = 2M\M 2 [uj ± 1) which 
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gives p = M2V '^ 2 — 1 where p is the momentum of the recoiling daughter baryon. The 
relevant expansion parameter close to the zero recoil point is thus y/u — 1. For example, 
at the zero recoil point only the helicity amplitudes H\ (allowed Fermi transition) and 

_ 2* 

H A = —H^ /y/2 (allowed Gamov-Teller transition) survive. In the LS-coupling scheme 

2 U 2 1 

with LS-amplitudes Tls, these correspond to the S-wave transition amplitudes TY± and 
TYi, respectively, where the orbital angular momentum L is defined with respect to the 

U 2 

relative orbital motion of the baryon B2 and the W ^_ shell in the rest frame of the baryon 
B\. In the literature one can find very ingenious approximation formulae for various decay 
distributions and polarization observables which are based on a near zero recoil expansion 
[HI [in]. These are usually referred to as effective theories of semileptonic hyperon decays. In 
this paper we shall, however, not discuss zero recoil or near zero recoil approximations but 
we always retain the full structure of the physical observables without any approximations. 

In order to get a feeling about the size of the helicity amplitudes we make a simple 
minimal ansatz for the invariant amplitudes at zero momentum transfer using SU(3) sym- 
metry. The analysis is greatly simplified by the fact that the C.G. coefficients for the 
(n — > p)-transition are the same as those for the (H° — > S + )-transition. One thus has 
F^(0) = 1 and F^(0) = 1.267 where the vector form factor -^(O) is protected from first 
order symmetry breaking effects by the Ademollo-Gatto theorem [391BU] . For the magnetic 
form factor F 2 V (0) we take -^(0) = M s o(/i p — (j, n )/(2M p ) = 2.6 as in [7j, where \i v and \i n 
are the anomalous magnetic moments of the proton and neutron. The second class current 
contributions are set to zero, i.e. we take F 3 y (0) = F 2 A (0) = 0. Note that a first class 
quark current can in principle populate the second class form factors and when 
Mi 7^ M 2 . For example, in the covariant spectator quark model calculation of [SI |4"5] one 
finds F 3 V (0) = (Mi - M 2 )/(6M 2 ) = 0.0176 and F 2 A (0) = 0. However, since we are not 
including SU(3) symmetry breaking effects for the other form factors we set F^ and F% 
to zero for consistency reasons since they vanish in the SU(3) symmetry limit. For ^3(0) 
we use the Goldberger-Treiman relation F 3 A (0) = M H o(M H o + M s+ )F 1 A (0)/(m^-) 2 = 17.1 
(see e.g. [3T]). 

For the g 2 -dependence of the invariant form factors we take a Veneziano-type ansatz 
which has given a good description of the g 2 -dependence of the electromagnetic form factors 
of the neutron and proton [16] . We write 

n=0 1 ~ m^W-^ V n=O m V,A + Ua J 

For Fi (q 2 ) and F^iq 2 ) we use ray = m#-*(892) — 0.892 GeV which is the lowest lying 
strange vector meson with J p = 1~ quantum numbers. Correspondingly we use itla = 
m K * { i.27o) = 1.273 GeV (J p = 1+) for F A (g 2 ) and m A = m K = 0.494 GeV (J p = 0~) for 
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(m B o - m. s+ ) 2 
q 2 [GeV 2 ] 



Figure 1: The g 2 -dependence of the six independent helicity amplitudes S(V, A) 

t-^ y A ' i-^r A i-rr A j-rr 



\fH v L f,T(V,A) := yftpHXf and L(V, A) := y^ H If multiplied by y^q 2 for the e-mode 

2 2 2 U 

(full range) and for the /x-mode (to the right of the vertical line q 2 = m 2 ). We also display the 
corresponding helicity amplitudes with the contributions of F^iq 2 ) and F^(q 2 ) switched off (dashed 
lines). 



V.A 



F A (q 2 ), respectively. The slope of the Regge trajectory is taken as a 1 = 0.9 GeV" 2 . The 
number of poles in Eq. (fT4|) is determined by the large q 2 power counting laws. One thus 
has ni — 1 and n 2 3 = 2. For the slopes of the form factors we thus have 1.781, 2.113, 
0.983 and 5.241 GeV -2 for FY , FY , F A and F A , respectively. These values correspond to 
("charge") radii (r 2 ) of 0.416, 0.494, 0.230 and 1.224 fm 2 , respectively, for the four above 
form factors. The ("charge") radii of the FY and F A form factors are close to the radii 
calculated in [47] using chiral input and a relativistic constituent quark model. The in- 
dependence of the form factors introduces only small effects since the range of q 2 is so small 
for the S° — > S + transitions. For the largest g 2 -value = (Mi — M2) 2 , the form factors 
have only increased by 2.8% {FY), 3.3% (FY), 1.5% (F A ) and 8.2% (F.f) from their q 2 = 
values. 

Based on these estimates for the invariant form factors we show in Fig. [T] a plot of the 
g 2 -dependence of the six helicity amplitudes. For easier comparability we have plotted 
the quantities ^/q^ 2 H^^ w . Close to the lower boundary q 2 = m 2 the longitudinal and 
scalar helicity amplitudes dominate, with iff ~ HY. and H A ps H A . Close to the upper 

2 U 2 2 U 2 

boundary at the zero recoil point q 2 = (Mi — M2) 2 , which is the relevant region for the 
/i-mode, the orbital S*-wave contributions HY. and = —\/2H A are the dominant 

2 2 2 U 

amplitudes with HY = -(F± /F A )Hf when F 3 y = F A = 0. In Fig. CD we also show 

2* 2 U 

plots of the helicity amplitudes where the contributions of the form factors FY (q 2 ) and 
F A (q 2 ) have been switched off. The difference is not discernible at the scale of the plot 
for the helicity form factors HY n , HY., H a and H A and is small for H A . The difference 

2 U 2 2 2 U 2 

is, however, sizeable for HY,- This can be understood from Eq. ([9]) which shows that, 
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even though H\ is small because it is proportional to y/Q- and thus of overall order 

2 1 

0((Mi — M 2 )/Mi), the relative contribution of F% to HY, is not suppressed by a factor of 

2 1 

0((M 1 — M 2 )/Mi) as it is in the other helicity form factors. Thus a measurement of the 
helicity form factor HY, would be ideally suited to determine the strength of _F 2 V . In fact, 

2 1 

in Sec. S] we shall discuss a forward-backward asymmetry measure which is proportional 
to HY, and is thus well suited for a measurement of F% ' ■ 

2 1 

We caution the reader that our ansatz for the form factors is only meant to implement 
the gross features of the dynamics of the semileptonic hyperon decays 5° — > S + + l~ + v\ 
which will eventually be superseded by the results of a careful analysis of the decay data. 
We shall nevertheless use the above minimal model for the S° — > S + form factors to 
calculate branching rates, the rate ratio r(e)/r(/x) and a forward-backward asymmetry in 
Sec. HI the longitudinal and transverse polarizations of the daughter baryon and the lepton 
in Sec. [5] and a mean azimuthal correlation parameter in Sec. [6] for this decay. In Sec. [7] 
we also show some Monte Carlo plots which are again based on the minimal form factor 
model. 

In the subsequent sections the rate and the angular decay distributions will mostly be 
written in terms of bilinear products of the sum of the vector and axial vector helicity 
amplitudes 

H\ 2 \ w = H^ 2Xw + H^ Xw , (15) 

since it is this combination which appears naturally in the master formulas describing the 
rate and the various decay distributions. One can of course rewrite the rate and the decay 
distributions in terms of bilinear products of the vector and axial vector helicity amplitudes 
FtX 2 \ w and . For the decay distributions this can be quite illuminating if one wishes 

to identify the overall parity nature of the observables that multiply the angular terms in 
the angular decay distributions. 



3 Unpolarized decay rate 

The differential decay rate is given by (see e.g. [2 



dT _ G\ 2 1 
dqHE l (2tt)3 |Kms| %M\ ^ ' 1 ] 

where is the usual lepton tensor (£0123 = +1) 

2 2 

IT = p? P : + rftf - -^9^ ± ie^ VhaVv3 . (17) 

As stated in the introduction the upper sign refers to the v{) case, whereas the lower 
sign refers to the v{) case. 
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The hadron tensor is given by the tensor product of the vector and axial vector 
matrix elements defined in Eqs. (CD) and (j2J), cf. 

= (M v + M A ), (M v + M A )l . (18) 

Eq. ([16]) shows that L^H^ determines the dynamical weight function in the (q 2 , E{) Dalitz 
plot (see e.g. the discussion in [43J ) . In a Monte Carlo generator one would thus have to 
generate events according to the weight L^ U H^ U in the (q 2 , E{) Dalitz plot. 

The differential g 2 -distribution can be obtained from Eq. (fl6l) by ^-integration, where 
the limits are given by (see e.g. [2] 



Et = ^2 (W + m 2 ) ± p(q 2 - m 2 )) . (19) 

Finally, in order to get the total rate one has to integrate over q 2 in the limits mf 
q 2 <: (M 1 - M 2 ) 2 . 

On reversing the order of integrations, the differential lepton energy distribution can be 
obtained from Eq. ( fl6l) by g 2 -integration. The relevant integration limits can be obtained 
from the inverse of Eq. ffT9]) . One obtains (see e.g. |24j ) 

q% = -(b±Vb 2 -ac) , (20) 

(X 



where 



M 2 + mf - 2 Mi Ei , 

M^M 2 - M 2 + m 2 - 2M 1 E l ) + m 2 M 



2 ■ 



m 2 ((M 2 - M 2 ) 2 + m 2 M 2 - (M 2 - M 2 )2M 1 E l 



Using Eiimax) := E milx = (M 2 — Mf + m 2 )/(2Mx) and _E;(min) := E min = mi one can 
simplify Eq. ( |20l to write 



ql = 2M 1 (E m l M - 1 Ei) + M 2 ( (Emax " El) ( El ± \/ Ef ~ E - ) + W) • (21) 

Finally, in order to get the total rate, one has to integrate over the lepton energy in the 
limits mi^Ei^ (Ml - Mf + m 2 )/2M 1 . 

As it turns out the two-dimensional integration becomes much simpler if one considers 
the two-fold differential rate with respect to the variables q 2 and cos 9 instead, where 9 is 
the polar angle of the lepton in the (I, vi) cm. system relative to the momentum direction 
of the Woff-sheii- Ei and cos# are related by (see e.g. 



2g 2 ^-g (g 2 + m 2 ) 

cos 9 = — ^ • ( 22 ) 

p(q z — mf) 
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Differentiating Eq. fl22l) one has 

d cos 9 2q 2 



dEi p(q 2 — m 
which leads to the differential decay distribution 



(23) 



dqH cos 9 ~ (2tt) 3 1 us| WM 2 q 2 ' 1 J 

It is clear from comparing Eqs. ffTBl) and (liHl) that, when writing a Monte Carlo program, 
one should not generate events in the (q 2 , cos#) Dalitz plot according to the weight L^ V H^ U . 

The cos 9 dependence of L^ V H^ U can be easily worked out by following the methods 
described in |24j which is based on the completeness relation for the polarization four- 
vectors 

£ ^Hf w K)w = f. (25) 

m,m'=t,±,0 

The tensor g mm i = diag(+, — , — , — ) is the spherical representation of the metric tensor 
where the components are ordered in the sequence m, m! — t, ±1, 0. One can then rewrite 
the contraction of the lepton and hadron tensors L^H^ as 



m,m' ,n,n' 




We shall refer to the second and third lines of Eq. fl2T)|) as the semi-covariant representation 
of the angular decay distribution. 

One has to remember that Eq. (126]) refers to the differential rate of the decay of an 
unpolarized parent hyperon into a daughter baryon whose spin is not observed. This means 
that one has to take into account the additional conditions Ai = X'i and A2 = A' 2 in Eq. (1261) . 
Angular momentum conservation then implies that not all index pairs m — w! and n = n' 
in Eq. (|26|) can be realized. Taking angular momentum conservation into account one has 
diagonal contributions m = ml = n = n' = t, ±1, as well as nondiagonal contributions 
with m — m' — t and n = n! = and vice versa. 

The point of writing L tlu H flu in the factorized form of Eq. (1261) is that each of the two 
factors in the second line of Eq. (1261) is Lorentz invariant and can thus be evaluated in 
different Lorentz frames. The leptonic part will be evaluated in the (/, v{) center-of-mass 
(cm.) frame (or PVoff-sheii^rest frame) bringing in the decay angle 9, whereas the hadronic 
part will be evaluated in the S~ rest frame bringing in the helicity amplitudes defined in 
Sec.[2J 
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Turning to the (/, v\) cm system the lepton momenta in the (x, y, z)-system read (see 

Fig. ED 



pf = (E t ;pi sin6, 0,- Pl cos6), (27) 
Pu = Pi (i, — sin9,0,cos9) . 

The angle 9 is always measured with respect to the direction of the lepton I, regardless 
of whether we are dealing with the (l~,Pi) or the (Z + ,z^) case. Since the orientation in 
the (x, ?/)-plane need not be specified in the present problem we have chosen the lepton 
momenta to lie in the (x, z)-plane. E\ and pi are the energy and the magnitude of the 
three-momentum of the charged lepton in the (I, vj) cm system, respectively, given by 
Ei = {q 2 + m 2 )/2^/lf and pi = (q 2 — mf) /2y/(p ■ The longitudinal and time-component 
polarization four-vectors take the form e M (0) = (0; 0, 0, —1) and e M (t) = (1; 0, 0, 0) whereas 
the transverse parts are unchanged from Eq. (jSJ). Using the explicit form of the lepton 
tensor Eq. (1171) it is then not difficult to evaluate Eq. (12 6j) in terms of the helicity amplitudes 
H\ 2 \ w of Sec. El One obtains 



o 

-(q -m 



-(1 =F cos^) 2 !^^! 2 + -(1 ± cos6') 2 |ii"_i_ 1 | 



+ i sm 
^ 2 '3 



\H 



\H 



+ | cos 2 9{\H ¥ \ 2 + \H_ ¥ \ 2 ) - 3 cos9(H ¥ H, + H_i t H_i Q ) 



(28) 



where the H\ 2 \ w are the sums of the corresponding vector and axial vector helicity am- 
plitudes defined in Eq. ( 1151) . We mention that the helicity flip factor mi/2q 2 does not give 
rise to a singularity since q 2 > mf. 

By explicit verification, or by hindsight, one can show that Eq. (T26|) can be written very 
compactly in terms of Wigner's d J -f unctions. One has what we shall refer to as our first 
master formula 



1 



E ( 

Ai,A2,Ayi/,A^ / ,J,J' 



-1) 



J+J'\U 



A; A„— I 



'X 



x 5 



• di 



A 2 Avi/-^A 2 A' 



(29) 



Except for the phase factor (— 1) J+J ' the master formula can in fact be derived by repeated 
application of the basic two-body decay formula in Appendix |A] The Kronecker 5-function 
S\ a -\ W) x 3 -\' in Eq. (1291) expresses the fact that we are dealing with the decay of an 
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unpolarized parent hyperon. One has to remember that Xw = and Xw = t both refer 
to the helicity projection (see Sec. [2]). Therefore there are nondiagonal interference 
contributions between J = 1, Xw = and J = 0, Xw = t because they are allowed by the 
angular momentum conservation condition A 2 — Xw = A 2 — X' w implying Xw = X' w . The 
interference contributions carry an extra minus sign as can be seen from the phase factor 
(— 1) J+J ' in Eq. ( |29|) . The phase factor (— 1) J+J ' comes in because of the pseudo-Euclidean 
nature of the spherical metric tensor g mm i defined after Eq. (1251) . 

The sign change in the first line of Eq. (1281) going from the (l~ , v\) to the (l + , vi) case 
can now be seen to result from the products of the relevant elements of the Wigner's d 1 - 
functions. For example, for A 2 = 1/2, Xw = 1 the nonflip contributions (A; = — X u = =Fl/2) 
are proportional to (d\-^ 1 ) 2 = (|(1 =F cos#)) 2 . There are no corresponding sign changes in 
the other lines of Eq. (l28i) . 

The h\ l \ v are the helicity amplitudes of the final lepton pair in the (/, v) cm. system. 
For example, for the (l~,i>) case with pi- along the positive z-axis, they can be worked 
out by using Eq. (jlj), the negative energy spinor of the massless antineutrino with helicity 
A p = | given by 

v*(l) = VK( \ ) , (30) 
and the SM form of the lepton current (Xw = X t - — Xp) 

V. Ti ^^,-H)7"(l + 7>^){ e ^-; ( ) 0) }. (31) 

We shall refer to the upper case A;- = — | as the nonflip transition and to the lower case 
A;- = | as the flip transition. Note the unconventional form of the SM lepton current 
which is due to the 75 definition in Sec. [2j The polarization four-vectors are given by 
£ M(t) = (1; 0, 0, 0), 6^(0) = (0; 0, 0, 1) and e**(±l) = (0; ^1, -i, Q)/y/2. The flip contribution 
is identical for Xw = t and A^ = 0. A similar expression can be written down for the case 
(7 + , V]) which we shall not work out in explicit form. For the moduli squared of the helicity 
amplitudes one finally obtains 

nonflip (A w = Tl) : K=t§,a,=± \ f = " m i) > ( 32 ) 

9 

flip (A w = t, 0) : |h Al=± i iA „ =± i | 2 = 8^(q 2 - m]) . (33) 

The notation in Eqs. (!32l and (|33|) (and elsewhere in the paper) is such that the upper and 
lower signs refer to the configurations (l~ , v\) and (l + , ui), respectively. In Eq. (1291) the sum 
over J, J' runs over and 1 and the index Xw, X' w runs over the four components t, ±1, 0. 
As remarked on before one has to remember to include the interference contribution from 
(J = 0; Xw = t) and (J = 1; Xw = 0) giving an extra minus sign. The matrix rf^ m ,, finally, 
is Wigner's <i 1 -function (d° mm , = 1 for m, m' = t) listed in Appendix B. 
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The form Eq. (|29|) readily affords a physical interpretation. H\ 2 \ W H^ X , determines 
the density matrix of the Woff-sheii (which happens to be block- diagonal in the present 
application). The density matrix is then "rotated" into the direction of the lepton in the 
(/, v\) cm. system with the help of the c^-functions whence the squared helicity amplitudes 
|^a ; a„| 2 determine the helicity dependent rates into the lepton pair. 

Performing the sum in Eq. ([29]) (Aj = ±1/2; X w = t, ±1, 0; J = 0, 1; J' = 0,1; A 2 = 
±1/2) one recalculates Eq. (|28|) . Note that the flip contribution proportional to m 2 /2q 2 
and nonfrip contributions are clearly separated in Eq. (T28]) . This separation facilitates the 
determination of the longitudinal polarization of the lepton to be discussed in Sec. El 

The differential rate dT /dq 2 is obtained from Eqs. (!24l) and ( J28l) by cos ^-integration 
which, in a sense, is trivial. One obtains 



dq 2 3 (2tt) 3 1 usl 8M 2 q 2 



\h h \ 2 + iiy_i_.il 2 + |#i | 2 + |#_i | 2 

2 2 2 2 



2 

+ ^^(|^| 2 + |ii_^| 2 ) + |ii, 1 | 2 + |ii_,_ 1 | 2 4 \H^ + \H 



2 2" J 2 2 2 2 u 



(34) 



The remaining g 2 -integration (mf ^ q 2 ^ (Mi — M 2 ) 2 ) has to be done numerically because 
of the nontrivial g 2 -dependence of the invariant form factors. 

A check on Eq. (134"]) is afforded by recalculating the Standard Model (SM) formula for 
semileptonic free quark decay q x — > q 2 + / + v setting F^(q 2 ) = Ff(q 2 ) = 1 and F 2 V '^ = 
in Eq. (131]) . One obtains 

' /] " M ' T ^^4p( - 2? + q 2 (l + p 2 ) + (1 - p 2 ) 2 



dq 2 q 

+ ^{ - 2g 4 - 2g 2 (l + p 2 ) + 4(1 - p 2 ) 2 }) , (35) 

where we have introduced scaled variables according to p = p/Mi, q 2 = q 2 /Mf, p 2 = 
M 2 /M 2 and r/ 2 = mf/M 2 , and where p = ±(1 + p 4 + q 4 - 2p 2 - 2q 2 - 2p 2 q 2 ) 1 ! 2 is the scaled 
magnitude of the daughter baryon's three momentum in the rest frame of the parent baryon. 
Also we have introduced the Born term rate 

n G 2 F \V us \ 2 Mf 



1927T 



?> 



(36) 



which represents the Standard Model decay of a massive parent fermion into three massless 
fermions, i.e. Mi ^ and M 2 ,m/,m^ = 0. The result Eq. fl35l) agrees with the SM result 
given e.g. in [42] . 

The g 2 -integration of Eq. (135j) can now be done analytically. One obtains (rj 2 ^ q 2 ^ 
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pf) 



YSM 



-24p 4 (l-r ? 4 )ln( 1 + p2 ~ r/2 ~ L ) + (p <- „) 



2p , ■ , - f/J , (37) 

where L = (1 + p 4 + r? 4 — 2p 2 — 27/ 2 — 2p 2 r] 2 ) 1 ^ 2 . Eq. (1371) again agrees with the result given 
in e.g. [12]. The symmetrization p rj in Eq. (13~T|) must be done for both the logarithmic 
and nonlogarithmic terms. The symmetry of the rate expression in Eq. (1371) under the 
exchange (p <->• 77) reflects the simple Fierz property of the SM (V — A) coupling. We 
mention that a less symmetric form of Eq. ( 1371) has been written down in [4~8] . 

We conclude this section with a comment on the relative merits of the two equivalent 
decay formulas in Eqs. (126]) and (1291) . In the semi-covariant representation Eq. (126]) the 
origin of the phase factor (— 1) J+J ' is clearly identified. Also, Eq. ( 1261) does not depend on 
the phase conventions chosen for the polarization four-vectors since they always enter in 
squared form. This is different in the master formula in Eq. (129]) and the master formulas 
written down in the following sections. They depend on the correct choice of phases for the 
polarization four-vectors and for the matrix elements of Wigner's <i J -functions. Judging 
from the fact that there exist different conventions for these phases in the literature the 
reader can appreciate what a hazardous enterprise it can be to get all the signs correct 
in the angular decay distributions if one has to rely solely on master formulas without 
explication of phase conventions. Whereas the signs of the polar correlations can usually 
be checked by angular momentum considerations there is no easy way to check on the signs 
of the azimuthal correlations to be discussed in the subsequent sections. In fact, we have 
repeatedly used the semi-covariant representation Eq. ( 1261) to check on the correctness of 
the phase conventions for the polarization four-vectors and Wigner's <i J -functions used in 
the different master formulas in this paper. 

4 Some numerical results 
4.1 Total semileptonic rates 

In order to obtain the total semileptonic rate we numerically integrate the differential rate 
Eq. ( 1341) over q 2 in the range m 2 ^ q 2 ^ (Mi — M2) 2 using the minimal model form factors 
described in Sec. |2j For the e-mode one obtains 

-> £+ + e~ + v e ) = \V US \ 2 1.169 • 1(T 17 GeV (38) 

This translates into a branching ratio of 

Si^S -> S + + e~ + v e ) = 2.513 • 10" 4 (39) 
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where we have used central PDG values for the lifetime of the S° (t(H°) = 2.90 • 10~ 10 s) 
and the CKM matrix element V us (V us = 0.2255) |49j . The branching ratio ( 1391) agrees very 
well with the experimental branching ratio BR(S° — ► S + + e~ + z/ e ) = (2.53 ± 0.08) ■ 10~ 4 

For the p-mode we obtain 

T(S° -> £+ + p~ + z/ M ) = |Ks| 2 9.893 • 10- 20 GeV , (40) 

which corresponds to a branching ratio of 

5^(2° ^ £+ + p~ + z/ M ) = 2.127- 10~ 6 . (41) 

In the PDG listings one finds BR(~° -> £+ + ^" + z^) = (4.6 • 10^ 6 [49J. The 
NA48 collaboration has published a preliminary result on this branching ratio which reads 
(2.2 ± 0.3 ± 0.2) • 10~ 6 (102 events) [501 EE] using a much larger data sample than that 
of the KTeV Collaboration (nine events) [52] on which the PDG value is based on. Our 
branching ratio ( l4Tj) nicely agrees with the NA48 value but disagrees with the KTeV value 
by more than one standard deviation. 

As it turns out one can get quite close to the exact results on the total semileptonic rates 
in a simplified setting. The resulting formulas are quite useful for a first assessment of the 
dominant dynamics of semileptonic hyperon decays. First, one neglects the contributions of 
the form factors F%'^ since they contribute at most at 0((M\— M2) / M\) . For the remaining 
two form factors and we neglect the g 2 -dependence since their g 2 -variation in the 
allowed g 2 -range is is only of 0(2%) (see Sec. [2]). With these simplifying assumptions the 
g 2 -integration of the differential rate in Eq. ( 1341) can be done analytically. The result is 
again written in terms of the scaled variables p 2 = M\jM\ and if = mf/M 2 , and the form 
L = (1 + p 4 + r] 4 - 2p 2 - 2 V 2 - 2p 2 T] 2 y/ 2 . One has 

r/r = ~ (VxT^p, v) + \f?\ 2 r{-p, tj)) , (42) 

where 



R(p,v) = L{A + A{p^T 1 )- 

-24p 3 ((l - rj 2 ) 2 + p 2 + p(l - r? 4 )) \nE 

-24?7 4 (1 - p 2 )(l - p + p 2 ) ln£(p <-> 77) , (43) 



and where we have used the abbreviations 



A = ±-7p 2 -7p* + pS + 6 v 2 p 2 -7 v 2 p\ 
B = 2p(l - hrj 2 - 2r] 4 + 10p 2 - 5p 2 r^ 2 + p 4 ) , 

E = 1 + "\-f- L . (44) 
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Upon setting = F^ = 1 in Eq. (142]) one reproduces the SM rate in Eq. ([3j 

Numerically, one has BR(E° -> £+ + e~ + u e ) = 2.451 ■ 10~ 4 and BR(~° -> £+ + p~ + 

z^) = 2.061 ■ 10~ 6 using the approximate result Eq. (j4"2l which is off the exact result in 

Eqs. ( ]39l and (HTT) by —2.4% and —3.1%, respectively. These deviations from the exact 

model figures provide a measure of the importance of the form factors F)f and F^, and 

the g 2 -dependence of the form factors. 

Considering the fact that (M1—M2) is small compared to (M1+M2) we expand Eq. (142]) 

in the variable 5={M X - M 2 )/{M l + M 2 ) = (1 - p)/(l + p) where 5 = 0.0501 in the case 

discussed in this paper. One obtains (see also |54j) 



r/r c 



512 1 



^I 2 (l + 3^)r(x) 5 5 + 0(5 7 ) 



5 (1 + 5) 

= l(l-p)Hl + P) 3 \Fr\ 2 (^+^-^)r(x)+0((l-p) 7 ). (45) 

The coefficient of the remaining <5 7 -term and the higher order terms in the small 5 
expansion Eq. (j4"5j) do not have as simple a structure as the remaining terms in Eq. (j4"5j) 
but can be calculated to be quite small (0.1% and 0.3%, respectively, in the e- and //-mode 
in the present case). It is curious that the next-to-next correction to Eq. (14"5]) is 0(5 9 ). 
The function r(x) in Eq. fj4"5l) is given by 



VT 3 ^, .2 „ A^ 15 4l 1-VT 



r( x ) = v x - (2 - 9a; 2 - 8s 4 ) - ya; 4 In ^ — — , (46) 

where x = 77/(1 — p). The overall factor (1 + 5)~ 8 in Eq. (|45j) can be seen to arise from 
power counting in Eq. (1431) considering the fact that p and 77 are proportional to (1 + ( s >)~ 1 . 

The small 5-expansion Eq. (j4"5l) is quite remarkable on two accounts. First, it is quite 
accurate since the corrections to Eq. (145]) set in only at 0(5 7 ) and not at 0(5 6 ) as one 
would naively expect, i.e. the corrections to Eq. ( ]4"5l) are only of 0(5 2 = 0.00251). From 
Eq. ( 1431) one can quickly appreciate that the rate measurement is quite sensitive to the 
ratio F^/FY ■ Second, it is quite remarkable that Eq. (145]) factorizes into a lepton mass 
independent and lepton mass dependent term. This will be very useful for a quick estimate 
of the rate ratio r(e)/r(/i) to be discussed in the next subsection. 

4.2 The rate ratio r(e)/r(/i) 

Of interest is the rate ratio T{e)/Y(p) which has been measured by the KTeV collaboration. 
Dividing the e-mode rate in Eq. (138]) by the /i-mode rate in Eq. (1401 we obtain 

r(e)/r(Ai) = 118.13 (55.6±?I;|[52]). (47) 

We have added the corresponding published experimental rate ratio and its errors from 
the KTeV collaboration in brackets. The KTeV value is off by more than two standard 
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deviations from the calculated value in Eq. (1471) . We mention that the NA48 Collaboration 
cites a preliminary value of 114.1 ± 19.4 for the rate ratio T(e)/T(n) [5()1 EI] which is in 
very good agreement with the calculated result in Eq. (j4"Tj) . 

Using the simplified form in Eq. (1421) we calculate T(e)/T({x) = 118.96 which is quite 
close to the full model value in Eq. (1471) . 

Finally, we consider the small ^-approximation of Eq. ( )42l) given by Eq. ( )45l) . Since 
the small ^-approximation factors into a lepton mass independent part and a lepton mass 
dependent part r(x) the rate ratio r(e)/r(^i) is simply given by the ratio r(x e )/r(x tl ). In 
particular this shows that, at this level of approximation, the rate ratio does not depend 
on the actual values of ^^(O) and F^(0). This in turn implies that a measurement of the 
rate ratio T(e)/T(/j,) does not reveal much of the dynamics of semileptonic hyperon decays 
except for a test of e-/i-universality. Numerically, one obtains 

r(e)/I» = r(x e )/r(x^) = 118.71 , (48) 

which is quite close to the corresponding ratio of 118.96 calculated in the simplified set- 
ting according to formula Eq. ( l4"2l) as one would expect from the quality of the small 
^-approximation Eq. (|45|) . 



4.3 Forward-backward asymmetry 



A very useful measure is the forward-backward asymmetry of the lepton in the W fj_ s h c ii 
rest frame (or (/, v{) cm. frame) defined by 



A FB (q 2 ) 



dT/dq 2 (forward) — dT / dq 2 (backward) 
dT /dq 2 (forward) + dT / dq 2 (backward) 

= N(q 2 ) 
■ D(q 2 )' 

The numerator factor can be calculated from Eqs. (1241) and (|28|) and reads 



(49) 



N(q z 



(2nY 



K. 



mf) p 



8M 2 q 



2„2 



mf 



-rljV TjA O t ( uV ttX 



rV 



Hi t Hi ) 

2 1 2 U 



i.e. 



D(q 2 



The denominator factor is simply given by the differential rate Eq. 
dT/dq 2 . 

For mi = (which is realized for the e-mode for all practical purposes) the forward- 
backward asymmetry can be seen to be directly proportional to the helicity amplitude HY, 

2 

and is thus very sensitive to the ratio of form factors F% jFy as discussed in Sec. EJ In the 
/i-mode there is in addition some sensitivity to the helicity amplitude which implies 

2* 

a certain sensitivity to the form factor ratio F^/F^. Looking at the size and signs of the 
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helicity amplitudes in Fig. [T] one finds that the forward-backward asymmetry is positive 
for the S° — ► S + transition in the e-mode and negative for the /i-mode. 

When averaging Afb(q 2 ) over q 2 the q 2 integration has to be done separately in the 
numerator and the denominator of Eq. (14"9]) . Using again the minimal form factor model 
of Sec. |2] we obtain 

(A FB ) (e-mode) = 0.081 (0.014) , (50) 

where we have added the corresponding value with F% switched off in brackets. The big 
difference between the two numbers in Eq. ( 150]) underscores the sensitivity of the forward- 
backward measure to the form factor ratio F% j 'F\ . For the /i-mode one obtains 

(A FB ) (/i-mode) = -0.082 (-0.115; -0.087) , (51) 

where the two figures in brackets refer to F^ and i*g , respectively, being switched off. 
There still is a sensitivity to the ratio F^ /F^, but at a much reduced level compared to 
the e-mode. The second number in the brackets of Eq. (I5TT) indicates that there is a very 
small sensitivity to the ratio F^/F^. 



5 Single spin polarization effects 
5.1 Polarization of the daughter baryon 

The lepton-hadron contraction L^H^ given in Eqs. (1281) and (1291) can be separated into 



contributions of positive and negative helicities of the daughter baryon denoted by L^ U H ±± 
They are given by 



-(IT cos 9) 2 \H i x | 2 + - sin 2 0\H H 



+ Hi {*\Hi t \ 2 + -\Hi 1 \ 2 sm 2 e 
2q 2 \2 l 5* 1 4 1 2 11 



+ -cos 2 6\Hi \ 2 -3cos6Hi t Hi 

£ ^ 2 2 



(52) 



(q 2 - m 2 



'I ± cos 9) 2 \H _i,| 2 + -sin 2 0|F_i n | 

V / I 2 1 1 4 2 



mt 



+ — ^ < -\H iJ 2 + -\H i -.rsin 2 ^ 
2q 2 1 2 1 ^ 4 1 -2- 1 



3 
2 



-r ^cos 2 fl|#_i | 2 - 3cos6H_i t H_ lQ 



(53) 



This allows one to compute the component P z of the polarization vector along the 
direction of P2 in the rest system of B\. One obtains 



Pm{B) 



(54) 
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On averagcl one has (P z ) = —0.65 and (P z ) = —0.33, respectively, for the e~-mode and 
//"-modes. 

In a similar vein the polarization of the daughter baryon in the x-direction can be 
obtained from Eq. (1291) by leaving the helicity label A2 unsummed. The product of helicity 
amplitudes now reads H\ 2 \ W H^ X , and the Kronecker 6 turns into 5a 2 -a w , x' 2 -x' w because, 
again, the parent baryon is taken to be unpolarized. As before, Xw = t and Xw — have 
zero helicity but transform as J = 1 and J = 0, respectively. One obtains 



2L H^ v 

PM = L, u Ht{+ L ~ iff ' (55) 



where 



2L» u H?_{e)=- 2 -{q*-m\ 



3 

sin0(±l - cos6)Hi 1 H_i 



2y/2 



3 

H =sm9(±l + cos9)Hi n H 1 , 

2V2 5 ° ~ 



2 



(56) 



+ 2^{^ si ^ cos ^|i^o-^-|-i) 

Of course, if one does not define a transverse reference direction the specification of P x 
does not make physical sense per se. Such a transverse reference direction is e.g. provided 
by the transverse momentum of the lepton in the semileptonic decay. In fact, we shall see 
in Sec. [6] how the density matrix of the daughter baryon enters the joint angular decay 
distribution of the cascade decay S° — > S + (— > p + tt + ) + l~ + v\ where the transverse 
reference direction is defined by the decay S + — ► p + 7r + . The polarization component P y is 
zero because we assume that the invariant amplitudes and thereby the helicity amplitudes 
are relatively real. On average one has (P x ) = —0.57 and (P x ) = —0.17, respectively, for 
the e~-mode and //"-modes. 

5.2 Polarization of the lepton 

The lepton-side flip- and nonflip-contributions to L^ V H^ U are clearly identifiable as can 
be seen by an inspection of Eqs. ( 1291) and One can thus directly write down the 

longitudinal polarization of the lepton for the decay of an unpolarized parent hyperon at 
no extra cost. One has 

L^H^iflip) + L^H^(nonflip) ' 



2 When averaging P z over q 2 one has to separately integrate the numerator and denominator of Eq. (|54|) 
after restoring the factor p(q 2 — mf)/q 2 in both the numerator and the denominator. 
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For the decay S° — > £ + +Z + z// the longitudinal polarization of the electron is ~ —100% 
over most of the range of q 2 because m e w 0. This changes only for g 2 -values very close to 



the threshold q 1 



mi 



For the \x -mode the longitudinal polarization is quite small and 



negative and remains smaller than ~ 30% in magnitude over the whole g 2 -range as Fig. [2] 
shows. On average one has (Pi M ^) = —0.18. Judging from the fact that Pj' 1 ^ is small the 
helicity flip and nonflip contributions are of almost equal importance for the //"-mode. 

It is important to realize that the longitudinality of the polarization pjp is defined with 
respect to the momentum direction of the lepton in the (I, v{) cm. system and not with 
respect to the momentum direction of the lepton in the rest system of the parent baryon 
5°. If one needs to avail of the longitudinal polarization in the latter frame this can also 
be done using the helicity method as has been shown in [27J. 

As before, the transverse polarization of the lepton can also be obtained from Eq. ( f29l) 
by leaving the helicity label A; in Eq. (1291) unsummed. One then obtains the density matrix 
of the lepton which we write as (L^ lJ H tl ' y )x l \ ll ■ This allows one to extract also the transverse 
polarization of the lepton P^ ■ One obtains (see also [12]) 



2{L tlv H^). 



(58) 



In order to evaluate Eq. ( 1581) for the (l~ ,i>) case one needs the relation hii = 
\/ mf /2q 2 h_ii. In Fig. [2] we show the g 2 -dependence of the transverse polarization of the 

in the decay 2° — > + +Pi. The transverse polarization starts off at rather high pos- 
itive values close to g^ in = m 2 and drops to zero at the zero recoil point g^ ax = (Mi — M 2 ) 2 . 
For the e~ the transverse polarization is practically zero over the whole q ,2 -range. Because 
of the lack of structure in the e~-case we do not show a plot of the polarization of the 
electron. 
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Figure 2: Longitudinal and transverse polarization of the \i in the (fx 



cm. system. 
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Figure 3: Definition of the polar angles and Op, and the azimuthal angle \ describing the decay 
of a polarized H° using the lepton side as polarization analyzer. P denotes the polarization vector of 
the H°. The coordinate system (xi,yi,Zi) is obtained from the coordinate system (x,y,z) by a 180° 
rotation around the y-axis. 

5.3 Decay of a polarized parent baryon 

In this subsection we consider the decay of a polarized parent baryon and in turn determine 
the angular decay distributions of the leptonic side and the hadronic side relative to the 
polarization of the parent baryon. The polarization of the parent baryon is described by 
the density matrix 

1/1 + Pcos#p Psintfp \ 
P^K W ~ g ^ P S i n ^ p i _ P cos e P ) > ( 59 ) 

where we have assumed that the polarization vector of the parent baryon lies in the (x, z)- 
plane with positive x-component as shown in Figs. [3] and HI The rows and columns of the 
matrix in Eq. ( |59i) are labeled in the order (1/2, —1/2). 

5.3.1 Lepton side as polarization analyzer 

The angular decay distribution is a straightforward generalization of Eq. ( 1281) where one 
now has to include the density matrix of the decaying parent baryon B\. Also, the rotation 
of the density matrix of the Woff-sheii hrto the direction of the lepton now involves also the 
azimuthal angle x- This brings in the phase factor e i ( Aw '~ A w r )( 7r- x). The appropriate angle 
entering the phase factor is (tt — x) since the azimuthal angle has to be specified in the 
leptonic (xj, $/j)-plane (see Fig. [3]). Using Appendix IA1 one obtains the master formula 

W(9,x,9 P )cx PA 2 -A w ,A 2 -A W (^)(-l) J+J > A;A , =±1/2 re l(A -- A '- )( ^ ) x 

Ai,A2,Avi/,A^, J, J' 

x di w ^{0)d^_^{0)HM W H^ , (60) 

where \ v = ±1/2 (\ v = 1/2 for (l~ , Pi) and X u = -1/2 for (/+, v{)). 

Doing the helicity sums and putting in the correct normalization one obtains 

dq 2 d cos 6d X d cos 6 P 6(2vr) 4 ' ns| 8M?q 2 
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(61) 



A similar result was published in [33] . However, our result in Eq. flBTj) does not agree with 
the corresponding result in [33] . 

5.3.2 Hadron side as polarization analyzer 

Following the familiar procedure of building up the cascade decay in a quasi-factorized 
form one obtains the master formula 



2~\w,X 2 -X w \ U P) rL \2\w rL ^\ 



W(6 B ,<p B ,d P )<x (-l) J+J 'Px 

A;,Avi/,A^, J, J',A2,A2,A3 

x J\<f> t J 1 dcosB \ti XlK=±l/2 \ 2 e^~ x '^x 



H\ n \„ 7 Ht'y x 

2 A W 



x <, A; _ A ,(e)4;,A ; -A,(^)e l(A2 - A '^^I 2A3 (^)^A3(^) Kol 2 , (62) 

where the h B are the helicity amplitudes of the decay B 2 — > P3 + tt. Latter decay is as 
usual characterized by the asymmetry parameter 



a B 



2 U 


2 _ 


h B 

2 U 




2 + 


h B 1 



(63) 



The asymmetry parameter for the nonleptonic decay S + — > p + 7r° relevant to this paper is 
given by ag = — 0.980 ^0015 @S] - Note that the phase factor in Eq. (|62|) now is exp[i(A2 — 
A' 2 )0b] which is appropriate for the azimuthal angle (ftp measured relative to the (x, z)-plane 
(see Fig. gj). 
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Figure 4: Definition of the polar angles 9g and 9b and the azimuthal angle <pB in the cascade decay 
of a polarized rP using the hadron side as polarization analyzer. 

Doing the helicity sum and the integration in Eq. ( I62j) . and putting in the correct 
normalization one obtains 
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(64) 



, (65) 



where B(B 2 — > S 3 + 7r) is the branching fraction of the nonleptonic decay B 2 ^ B 3 + n. 



6 Joint angular decay distribution 

Following the familiar procedure the joint angular decay distribution for the semileptonic 
cascade decay B\ — > B 2 {— > B% + 7r) + £ + z/j of an unpolarized parent baryon Si can be 
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Figure 5: Definition of the polar angles 9 and Ob, and the azimuthal angle \ m the joint angular 
decay distribution of an unpolarized H° in the cascade decay H° — > p + 7r°) + l~ + P/. The 

coordinate system {x\, yi, z{) is obtained from the coordinate system (x, y, z) by a 180° rotation around 
the y-axis. 



derived from the master formul 

W(0,XM OC J2 (-1) J+J >U,= ± 1/ 2 I V(AW_AW)( ^ X) ^-A W -,A^A W X 

AjjAwjA^, J, J' ,\2,\ 2 ,\3 

X C^,-A.(f)<,^(e)fi-W^ 4 2 A 3 (^)4A3(^)l^ 3 0r • (66) 

The polar angles 9, 9b and the azimuthal angle x are defined in Fig. [5j The Kronecker 5 
in Eq. (166|) expresses the fact that we are dealing with the decay of an unpolarized parent 
hyperon which implies A 2 — Xw = X' 2 — X' w . 

When writing down the corresponding normalized decay distribution we shall as before 
assume that the helicity amplitudes are relatively real. One obtains 



dq 2 d cos 9d X d cos 9 B y A 7 6(2tt) 41 usl 8M 2 q 2 

as cos 9 b ) + -{1±cob6) 2 \H_i_ 1 \ 2 (1 — as cos9 B ) 
8 2 8 2 

+ - sin 2 9 (\Hi \ 2 (l + a B cos 9 B ) + \H_i Q \ 2 (l -a^ cos 9b) j 

± — i=aB sin 9 cos x sin 9 B f (1 T- cos9)H_i Hi 1 + (1 ± cos9)Hi Q H_ 

2y2 ^ 2 2 2 2/ 

m 2 f 3 3 
+ -^l-\Hi t \ 2 {l + a B cos9 B ) + -\H_i t \ 2 (l - a B cos9 B ) 

— 3 cos9yHi t Hi (l + as cos9b) + H_i t H_i (l — as cos#b) j 

+ -cos 2 6>(|iii | 2 (l + a B cos9 B ) + |#_i | 2 (l - a B cos9 B )^ (67) 



3 Much to the embarrassment of one of the present authors (JGK) there was a sign mistake in the 
azimuthal correlation term of the corresponding joint angular decay distribution written down in [28] for 
the semileptonic decay A c — > A(— > p + tt~) + l + + vi (mi = 0). The source of this error was that in 
[2"5] we used the phase factor exp [i(Aw — ^w)( — X)] instead of the correct form exp [i(Aw — A^ v )(7r — x)] 
to determine the sign of the azimuthal correlation term. This error was discovered and rectified through 
experimental evidence [55] . 
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(68) 



We have performed several checks on the correctness of the signs of the azimuthal 
correlation terms by using the semi-covariant representation Eq. (T2"6"l) and even doing a 
full-fledged covariant calculation^. The overall sign of the nonflip azimuthal correlation 
terms (sixth and seventh line in Eq. (l67j) ) corrects the sign mistake in [28]. Note the 
reciprocity of the angular decay distributions Eq. (IBTi) and Eq. ( 1671) . One obtains Eq. ( 1671) 
from Eq. ( IBTI) by the substitutions (1 + sign{\2 — \\v}P cos9p — > (1 + sign{\ 2 }aB cos9b) 
for the polar correlation terms and P sin 9pH\ 2 \ w H X2 \' w — > a B sm9BH\ 2 x w H_ x . 2 \> w in the 
azimuthal correlation terms. 

Eq. ( 1S"7j) can be cast into a form where the dependence on the polarization vector of 
the daughter baryon becomes explicit. One has 

^ -B(B 2 B 3 + *)\<%\ V ^*£g>LrH~x 



dq 2 d cos 9dxd cos 9 B v / 4(2vr) 41 usl 8M*q 2 

x (1 + P z a B cos 9 B + P x a B cos(7r — x) sin ^b) , (69) 

where L^H^, P z and P x are given in Eqs. (|28|) . (|54l) and (|55|) . respectively. When 
integrating Eq. (|69|) over cos#, cos9 B and g 2 one can define a mean azimuthal correlation 
parameter (7 ) through the relation T ~ 1 + (7 ) cosx- Using again the minimal form factor 
model of Sec. [2]one finds the numerical values (7 ) = —0.44 and (7 ) = —0.13 in the e- and 
/z-modes, respectively, for the mean azimuthal correlation parameter. 

At zero-recoil one finds a rather simple expression for the above azimuthal correlation 
parameter. It reads 

7zero recoil i „a " (^0 

!6 i + gL(i + 2|i7r ( |V|^ 1 | 2 ) 

which gives 7 = —0.61 and 7 = —0.17, respectively, in the e- and yU-mode in the minimal 
form factor model described in Sec. [2j Eq. (1701 shows that, in the e-mode and at zero 
recoil, the azimuthal correlation parameter is independent of the form factors as stated 
before in [28J. In the /x-mode, however, the azimuthal correlation parameter at zero recoil 
does depend on the form factors through the ratio HXjHf . Since HXJH^ = F 1 V '/ (\/2Ff) 



It is fair to say that the transcription of the results of a covariant calculation into the helicity frame 
results used in this paper takes considerable amount of calculational effort. 
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at zero recoil (assuming F, 



v 



F A 



0) this would afford the opportunity to determine 



the ratio FY jF\ through a zero recoil or near zero recoil measurement of the azimuthal 
correlation parameter in the /x-mode. 



7 Monte Carlo event generation and sample plots 

In this section we present a few sample distributions generated from our event generator in 
order to demonstrate the viability of our generator. As dynamical input for the form factors 
we used the minimal form factor model described at the end of Sec. [21 or slight variations on 
it. Of course, any other dynamical model can be used as input in the event generator. For 
the angular decay distribution we used the full five- fold decay distribution from Appendix [D] 
describing the full decay chain S°(|) — > > p + 7r°) + W~ s _ shell (— > F + Vj). Masses and 
the decay asymmetry parameter are taken from |49j. 
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Figure 6: Angle between the electron and the neutrino in the E° rest frame for different values of 
the ratio Ff(0)/FY(0). 
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Figure 7: Energy spectrum of the electron for different values of the ratio F^(Q)/FY(0) in the 3 
rest frame. 
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Figure 8: Left: Angle between S + and x axis in the H° rest frame. Right: Angle between proton 
and x axis in the S + rest frame for different initial state polarizations of the E° hyperon. 
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Figure 9: Comparison of lepton spectra for the e- and ^-mode 



The implementation was done as follows. We first generated the 3-body phase space 
of the primary decay H° — > S + + W oS _ shell using the widely used function genbod from the 
CERNLIB [72] library. Without loss of generality, the axis of the initial state polarization 
of the parent baryon was chosen to point along the lab x-axis. The momenta of the decay 
products of the secondary decays S + — > p + n° and W~ s _ shell — > l~ + \>i were generated 
with uniformly distributed directions. Since the secondary decays are two-body decays, the 
moduli of the respective momenta are fixed. The resulting momentum vectors were used 
to obtain the angles and momenta needed to calculate the value of the matrix element. 
This result was multiplied by the phase space factor of the primary decay returned by 
the genbod routine. Applying an acceptance-rejection method, the whole procedure was 
repeated until a generated event was no longer discarded. 

In Fig. [6] we show a plot of the dependence of the rate on the angle between the 
electron and the neutrino in the S° rest frame. In order to exhibit the sensitivity of this 
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distribution to the form factor ratio F^(0)/FY(0) we also show corresponding distributions 
with slightly varied form factor ratios. We mention that the distribution Fig. [6] and the 
following distributions are normalized to unity. 

In Fig. [7J we show the electron energy spectrum and its dependence on the form factor 
ratio F 2 y (0)/F 1 v (0). In Fig.[H]we show plots of the angular dependence of the angle between 
the E + and the x axis in the S° rest frame (left) and between the proton and the x axis in 
the S + rest frame (right) for different polarizations of the parent S° hyperon. In order to 
demonstrate the dependence on lepton mass effects, in Fig. [9] we show the energy spectra 
of the electron and the muon in the rest frame of the parent 2°. 

8 Summary and conclusions 

We have worked out the angular decay distributions that govern the semileptonic cascade 
decay H° — > S + (^ p + 7r°) + W^g_ she]1 (— > l~ + v\) using a cascade-type analysis. The 
cascade-type analysis has certain advantages, the main advantage being that one obtains 
the decay distributions in a compact quasi-factorized form. This leads to rather compact 
forms for the decay distributions. In our analysis we have included lepton mass effects as 
well as polarization effects of the decaying parent hyperon. We have always indicated the 
necessary sign changes when going from the (l~ , v\) case to the (/ + , v\) case. Our angular 
decay formulae are thus applicable also to the semileptonic hyperon decay S + — > A+e + +z/ e , 
or to semileptonic charm baryon decays induced by the transition c — > s + l + + v\ and also 
to the decays t — > b + Z + + v\ . It should be clear that our angular decay formula are 
also applicable to the corresponding nonleptonic baryon decays involving vector mesons 
(\w = ±, 0) or pseudoscalar mesons (A^ = t). In this case one has to omit the interference 
contributions between the time-component and the space-components of the currents. 

Of interest are also the corresponding semileptonic antihyperon decays. The angular 
decay distributions of semileptonic antihyperon decays can be obtained from the corre- 
sponding angular decay distributions of the semileptonic hyperon decays by the replace- 
ments H\ 2 \ W (B) — > H\ 2 \ W (B), cub — ► Oi B and changing from the to the (/ + ,z/^) 
case (or vice versa). Neglecting again CP-violating effects one has from CP-invar iance 
Hx 2) x w (B) = H_x 2 ,-\ w (-B) an d as — ~ a B- One can verify that the decay distributions 
in Eqs. ( l61i) . ( |64|) and ( 1671) are form invariant under H\ 2 \ w — > P_a 2 -a w , «b — ¥ —&b and 
P — > — P as follows from CP-invariance. We mention that the NA48 collaboration has 
recently observed the decay 5° — ► S+e + z/ e and, based on 555 events, have given a branching 
ratio of (2.55 ± 0.14 stat ± 0.10 syst ) • 1CT 4 for this semileptonic antihyperon decay [53] . 

We have summed over the helicity states of the final particles assuming that their 
polarization go unobserved. This corresponds to taking the trace of the density matrix of 
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the final particles. It is clear that one can equally well calculate the density matrix of the 
final state particles by leaving the relevant helicity index unsummed. This was illustrated 
for the density matrix of the final lepton in the semileptonic decay process. 

Doing the helicity sums in the master formulas listed in this paper by hand can become 
quite cumbersome. However, this task can be automated and can be left to the computer. 
The relevant Mathematica codes can be obtained from A. Kadeer. We mention also that 
the helicity frame analysis used in this paper can be easily transcribed to a transversality 
frame analysis (see e.g. [24J) where the z-axis is perpendicular to the hadron plane . In 
fact, any choice of z-axis in the analysis will provide the same total amount of information 
on the dynamics of the process entailed in the invariant amplitudes. It is then a question 
of experimental exigency of whether to analyze angular decay distributions in the helicity 
frame or the transversality frame, or, for that matter, in any other frame. 

For the sake of conciseness we have written our results in terms of bilinear products of 
the helicity amplitudes H\ 2 \ w = H XiX + H X2X instead of bilinear products of the vector 
and axial vector helicity amplitudes H X2X and H X2X . Writing the decay distributions in 
terms of H X2X and H X2X can be quite illuminating if one wishes to identify the overall 
parity nature of the observables that multiply the angular terms in the angular decay 
distributions. 

We have formulated a minimal form factor model assuming SU(3) symmetry at zero 
momentum transfer and a canonical g 2 -dependence of the form factors. We have used this 
minimal form factor model to numerically calculate a few observables such as branching 
rates, the rate ratio T(e)/T(n), a lepton-side forward-backward asymmetry, the longitu- 
dinal and transverse polarizations of the daughter baryon and the lepton and a mean 
azimuthal correlation parameter in the decay S° — > p + tt°) + l~ +ui. 

We have written a Monte Carlo event generator which is based on the angular decay 
distributions derived in this paper. Among others, the event generator allows one to gener- 
ate events in the parent baryon rest frame. The MC program can be obtained from Rainer 
Wanke of the NA48 Collaboration (Rainer.Wanke@uni-mainz.de). We have presented a 
few decay distributions and correllations based on this event generator. We have, however, 
not systematically investigated which observables would be optimal to obtain the maximal 
possible information on the underlying dynamics encapsuled in the invariant form factors 
or the helicity amplitudes. A discussion of these issues also using parent baryon rest frame 
decays can be found in e.g. [TJ. 

We did not provide an in-depth analysis of the dynamics of semileptonic hyperon de- 
cays as is necessary if one wants to extract a value of the CKM matrix element V us from 
semileptonic hyperon decay data. This issue was discussed in |56j E3 EH1 EH] . We mention 
that there has been some recent progress in the dynamical description of strangeness chang- 
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ing semileptonic hyperon transition form factors from the lattice community [601 ED E2] 
and in the framework of chiral symmetry [ITJ [66] including the explicit calculation of chiral 
corrections [63J, Ell ES] . 

We finally emphasize that we have not included CP- violating effects or radiative cor- 
rections in our analysis. The latter can be included using the results of [7j. It will be 
interesting to find out how the radiative corrections will affect the angular decay distribu- 
tions. 
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A Two-body decay of a polarized particle in the he- 
licity formalism 

In deriving the two-body decay of a polarized spin J particle in the helicity formalism we 
shall closely follow the approach of [HZl EE] ■ Consider the two particle decay a — > b + c of 
a spin J a particle where the polarization of particle a in the frame (xq,i/q,Zo) is given by 
p° XaX , . Consider a second frame (x, y, z) obtained from (x , yo, Zq) by the rotation R{6, <fi, 0) 
and whose z-axis is defined by particle b. The polarization density matrix p in the frame 
(x,y,z) is obtained by a "rotation" of the density matrix p° from the frame (xq, yo, z ) 
to the frame (x,y, z). The rate for a — > b + c is then given by the the sum of the decay 
probabilities | i^A 6 A c | 2 (with A a = \ — A c ) weighted by the diagonal terms of the density 
matrix p of particle a in the frame (x, y, z). Thus we find 

r a _ b+c (M)oc J2 \ H ^\ 2 K,x h -xS e ^)pl a ^AM-xS^^ ( A1 ) 

where 

D J m)m ,(e,<t>)=e- im *d J mm ,(0) . (A2) 
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All the master formulas written down in this paper can be obtained by a repeated 
application of the basic two-body formula Eq. (lAljl . 



B Wigner's d J — functions for J=l/2 and J=l 



For definiteness we list explicit forms for the two (i J -functions of Wigner used in this paper. 
We use the convention of Rose |69| which is also the convention of 149 1 . One has 
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(B2) 



for J = 1. The rows and columns are labeled in the order (1/2,-1/2) and (1,0,-1), 
respectively. 



C T— odd contributions 

In the main text we have assumed that the invariant form factors and thereby the helicity 
amplitudes are relatively real. If one allows for relative phases between the helicity am- 
plitudes one will obtain so-called T-odd contributions in the angular decay distributions. 
They appear in the azimuthal correlation terms as can be seen by the following example 
taken from the joint angular decay distribution in Sec. [61 One of the azimuthal correlation 
terms derives from the helicity configurations (A2 = —1/2, Xw = 0; A' 2 = 1/2, X' w = 1) and 
(A 2 = 1/2, X w = 1; X' 2 = — 1/2, A'^ = 0). Picking out the relevant terms in the master 
formula Eq. (|66|) one has 

ffiifT i0 e <(w-x) + H_i o m ie - i{n - x) = -2cosxRe(Hi 1 H% ) -2sinxIm(Hi 1 H% ) . 

22 (CI) 

The cosx dependent term already appears in Eq. (IBTl) whereas the sinx dependent term 
has been dropped in Eq. (IBTl) because of the relative reality assumption for the helic- 
ity amplitudes. Adding the relevant and 9b dependent trigonometric functions in 
the above azimuthal correlation term one has the two angle dependent T-odd terms 

(sin sinx sin 0# Im^i^H^ 1 )) and (cos sin sinx sin 9b ^(Hij^H^ 1 )) proportional to 

2 2 2 2 

sinx- 
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Next we rewrite the product of angular factors in terms of scalar and pseodoscalar 
products using the momentum representations in the (x, y, z)-system (see Fig. [5]). For the 
normalized momenta one has (p 2 = 1) 



Pi- = (sin# cosx, sin 9 sin x, — cos 9) , 

p w = (0,0,-1), 

p s+ = (0,0,1), 

p p = (sin ^,0, cos 6> s ) , 



(C2) 



where the momenta have unit length indicated by a hat notation. The above T-odd angular 
factors can then be rewritten as 



Under time reversal (t — > —i) one has (p — > —p). Since the T-odd momenta invariants in 
Eq. ( 1C3I) involve an odd number of momenta they change sign under time reversal. This 
has led to the notion of the so-called T-odd obsevables. Observables that multiply T-odd 
momentum invariants are called T-odd obsevables. They can be contributed to by true 
CP-violating effects or by final state interaction effects unless either or both change all 
helicity amplitudes by the same common phase. One may distinguish between the two 
sources of T-odd effects by comparing with the corresponding antihyperon decays since 
phases from CP- violating effects change sign whereas phases from final state interaction 
effects do not change sign when going from hyperon to antihyperon decays. 

From the above example it should be clear how to obtain the T-odd contributions 
from the master formulas for the other cases. In practise what one has to do is to add 
terms where the real part of the bilinear forms of helicity amplitudes is replaced by the 
corresponding imaginary part and the cosine of the azimuthal angle is replaced by the sine 
with a possible sign change. 

D Full five-fold angular decay distribution 

In this appendix we write down the full five-fold angular decay distribution for the semilep- 
tonic cascade decay of a polarized hyperon. There are now altogether three polar angles 
9, 9b and 9p, where 9p describes the polar orientation of the polarization vector of the 
parent hyperon as shown in Fig. [10] (which is directly taken from [IE]). Since there are 
now two planes in the cascade decay, there is one more azimuthal angle which we choose 
as (pi as shown in Fig. [TTJ It is important to note that Fig. [11] shows a special configuration 



sin 9 sin x sin 9b 
cos 9 sin 9 sin x sin 9 b 



Pw (Pi- x Pp) , 

(pi- ■ Pw) [Pw ■ (Pi- x Pp)] ■ 



(G3) 
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Figure 10: Definition of the three polar angles 9, 6b and 6p in the semileptonic decay of a polarized 
S° into E + + l~ +&i followed by the nonleptonic decay S + —>■ p + ir°. The polarization vector of the 
parent baryon P lies in the (x,z)-plane with positive P x component. 



where the momentum of the proton lies in the first quadrant and the momentum of the 
lepton lies in the second quadrant. It is clear that, for this special configuration, the three 
azimuthal angles 4>ii4>B and x a dd up to n (cpi + 4>b + X = n )- F° r other configurations 
it may happen that the three angles add up to tt + mod(27r) if the rotation sense of the 
angles in Fig. [11] is kept. This will be of no consequence for the angular decay distribution 
which is invariant under azimuthal 2tc shifts. 

The full five-fold angular decay distribution can be directly taken from [TB] after in- 
cluding the appropriate sign changes going from the (7 + , vj) to the vj) case§. We have 
simplified the corresponding expressions in [18] by assuming as before that the helicity 
amplitudes are real. For completeness we shall also write down the decay distribution in 
explicit form using Wigner's <i J -functions as before. One has the master formula 

w(o,e P ,o B ,<i>B,<in) oc Yl (-i) J+J >l Wl/2 |V (A -- A '^x 

x P\ 2 -\ w ^ 2 -\' w {6p)di w , Xl - Xi X6)di, w Xi _^{6)H X2Xw Hl, 2 x w X 

x e^-^^M^WlhU 2 ■ ( D1 ) 
For the normalized five-fold angular decay distribution one finds 

5 Apart from listing angular decay distributions Ref. [18] contains much additional useful material like 
e.g. a discussion of the statistical tensors of the processes and their bounds, HQET results for heavy 
baryon transition form factors, etc.. 
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Figure 11: Definition of the three azimuthal angles fa, 0b and x {'Pi+'Pb+X = 7r ) i n tne semileptonic 
decay of a polarized S°. Fig. [TT] is a view of Fig. [10]from the right along the negative z-direction. pf 
and p£ denote the transverse components of the momentum of the lepton and proton, respectively. 
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x 6™ + 3 cos 6gS + cos0 B C 



x 



(|coB a 6»-i) C 

— 3-\/2 sin cos 0; 6q} 

— 2 sin^ B cos 4>b &10 

3 

r= sin 26* cos 0; 6q 2 

v 2 

+ 3 cos 6* cos Ob &00 

+ (_COS 2 0--)COS0 B &00 

3 

— - \pl sin 9 sin # cos x b 
3 

— -v^sin^B sin 26* cos x b\\ 
3 

+ - sin 2 9 sin B cos(x - fa) b\l 

— 6 cos 9 sin 9 b cos0b b 10 

— 3\/2 sin 6* cos 6*b cos 0; b^\ 
3 

sin 29 cos # b cos fa 6q 2 

— sin 9 B {2> cos 2 6* — 1) cos0b&io 



(D2) 



It is important that the rotation sense of the azimuthal angles in Fig. [TT] is kept. We 
have used the relation fa + <p B + X — 71 + mod(27r) to rewrite cos(0b + 0/) = — cosx and 
cos(0b + 20;) = — cos(x — fa). Note that Eq. (ID2I) contains the redundant angle 0b- As 
before one can reexpress cos 0b as cos 0b = — cos(0z + x)- 
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The coefficients b^j in Eq. (1D2I) are given by § 

C=((l + ^)(l^-i-i| 2 + l^|or) + 3e|^ 
+ ((l + e )(|^ 1 | 2 + |i/„i | 2 ) + 3 e |iJ_ | 

C =l(T\H. 2l \ 2 -4eH_ h0 H 



Pa 

2 2 



2 1 / 2 2 



-i( T |^_ 1 | 2 + 4e^ 

C =a fl ((l + e)(-|^-|-x| 2 + l^ol 2 ) + 3e l^/)p|i 
- + e)(|^_ |0 | 2 -|^x 1 | 2 )+3e|^_x 

C =^(-2|F, | 2 + |^ 1 



P-i-i 

2 2 



b°{ =\(2eH_^ t H_ h _ 1 ±H_j l0 H. h _ 1 



b\° = - a B {?,eH ¥ H_ ¥ + (1 + e)Hi H_i ) p_ii , 
Kl =^(H_ h0 H_,_ 1 -H kl H^)p h _ k , 

^~(l-2e)(2|iJ_, | 2 + |^ 1 | 2 ) p_i_, 
-f(l-2 e )(|^_ 1 i 2 + 2|^ | 2 )^, 

b\\ =a B (2eH, l H_ ¥ TH ll H_ ¥ ) p_i_i_ 

+ a fl ( - 2etf |f #_ i _ 1 =F ^-i-i) P±i , 

6j 2 =aB(l-2e)(# §1 #_i p_i_i - Hi Q H_i_ lP i_i 
6i| = -a B (l-2 e )ffi 1J H-_i_ 1 pi_i, 

& io =«b e [Hi H_i t + Hi t H_i 



,11 _«B 

Ooi- 2 



(-2eH_ it ff_i_ lT ff_ Jo ff_j_ 1 



2tH ii H i. i H i-iH \ n ) Oi _i 



6 The coefficient 6} J takes twice the value as compared to the corresponding coefficient in [18]. Also in 
Eq. (48) of [18] concerning the overall normalization one has to effect the replacement g 2 — > (g 2 — mf) 2 /q 2 . 
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bH = - <f (1 - 2c) (ff-^-i-! + ^oj n-,- 



--a B (l-2e)Hi H_i p_ 



1 1 

2 2 



(D3) 



We have introduced the abbreviation e = mf /2q 2 for the leptonic flip suppression factor. 
As in the main text the upper signs in the coefficients b%j hold for the case (l~ , T>i) relevant 
to the cascade decay H° — ► p + 7r°) + l~ + v x treated in this paper. The lower signs 

hold for the case (/ + , v{) discussed in [H]. Finally, p Xl \' is t ne spin density matrix of the 
parent hyperon given in Eq. (1591) . 

We have performed various checks on Eq. (1D2I) . First we found it to agree with the 
angular decay distribution derived from the master formula Eq. ( 1D1I) . We further checked 
that Eq. (1D2I) reduces to the decay distributions listed in the main text after integration 
or after setting the relevant parameters to zero. We thus checked that Eq. (ID2I) reduces to 
Eq. (|67|) when setting P = 0. There is a factor of Air from the integration over cos 9p and 
4>i. Further Eq. ( 1D2I) reduces to Eq. (16TT) when setting olb = 0, dropping the branching 
ratio factor B(E>2 — > -B3 + M) and replacing (pi by (tt — x). Also there is a factor 47r from 
the integration over cos## and 0b- Finally, Eq. (1D2I) reduces to Eq. (16^1) when integrating 
over 0z and cos#. As mentioned before we have assumed that the helicity amplitudes (or 
the invariant amplitudes) are relatively real. Nonzero relative phases between the helicity 
amplitudes could arise from final state interaction effects or from extensions of the SM that 
bring in CP-violating phases (see e.g. [7Ql[7l]). In such a case one would have to keep the 
full phase structure contained in the master formula Eq. ( 1D1I) or in the original version of 
Eq. (|D2J listed in [33]. 
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